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^ ! Abstract 
' 

Q_i| One plus two-body embedded Gaussian orthogonal ensemble of random matrices with parity 

■ [EGOE(l+2)-7r] generated by a random two-body interaction (modeled by GOE in two particle 

spaces) in the presence of a mean-field, for spinless identical fermion systems, is defined, generalizing 
the two-body ensemble with parity analyzed by Papenbrock and Weidenmiiller [Phys. Rev. C 78, 
054305 (2008)], in terms of two mixing parameters and a gap between the positive (vr = -1-) and 
negative (vr = — ) parity single particle (sp) states. Numerical calculations are used to demonstrate, 
using realistic values of the mixing parameters appropriate for some nuclei, that the EG0E(l-|-2)- 

(N 

J> ' vr ensemble generates Gaussian form (with corrections) for fixed parity eigenvalue densities (i.e. 



state densities). The random matrix model also generates many features in parity ratios of state 
densities that are similar to those predicted by a method based on the Fermi-gas model for nuclei. 
We have also obtained, by applying the formulation due to Chang et al [Ann. Phys. (N.Y.) 66, 
137 (1971)], a simple formula for the spectral variances defined over fixed-(mi, m2) spaces, where 
mi is the number of fermions in the +ve parity sp states and m2 is the number of fermions in the 
— ve parity sp states. Similarly, using the binary correlation approximation, in the dilute limit, we 
have derived expressions for the lowest two shape parameters. The smoothed densities generated 
by the sum of fixed-(mi, 7712) Gaussians with lowest two shape corrections describe the numerical 
results in many situations. The model also generates preponderance of -|-ve parity ground states 
for small values of the mixing parameters and this is a feature seen in nuclear shell model results. 
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I. INTRODUCTION 



Random matrix theory (RMT), starting with Wigner and Dyson's Gaussian random 
ensembles [1, 2] introduced to describe neutron resonance data [3, 4], has emerged as a 
powerful statistical approach leading to paradigmatic models describing generic properties 
of complex systems [5-9]. Developments and applications of RMT in nuclear physics in 
last 30 years have been reviewed recently by Weidenmiiller and collaborators [10, 11]. The 
Wigner-Dyson classical Gaussian orthogonal (GOE), unitary (GUE) and symplectic (GSE) 
ensembles are ensembles of multi-body interactions while the nuclear interparticle interac- 
tions are essentially two-body in nature. This together with nuclear shell model examples 
led to the introduction of random matrix ensembles generated by two-body interactions in 
1970-1971 [12, 13]. These two-body ensembles are defined by representing the two-particle 
Hamiltonian by one of the classical ensembles and then the m (m > 2) particle if-matrix 
is generated by the Hilbert space geometry. Thus the random matrix ensemble in the two- 
particle spaces is embedded in the m-particle iJ-matrix and therefore these ensembles are 
gcncrically called embedded ensembles (EEs). Simplest of these ensembles is the embedded 
Gaussian orthogonal ensemble of random matrices generated by two-body interactions for 
spinless fermion (boson) systems, denoted by EG0E(2) [BEG0E(2); here 'B' stands for 
bosons]. In addition to the complexity generating two-body interaction, Hamiltonians for 
realistic systems such as nuclei consist of a mean-field one-body part. Then the appropriate 
random matrix ensembles arc EE(l+2). The spinless fermion/boson EGEs (orthogonal and 
unitary versions) have been explored in detail from 70's with a major revival from 1995 and 
it is now well understood that EGEs model many-body chaos or stochasticity exhibited by 
isolated finite interacting quantum systems [9, 14]. Besides the mean-field and the two-body 
character, realistic Hamiltonians also carry a variety of symmetries. In many applications 
of EGEs, generic properties of EGEs for spinless fermions are 'assumed' to extend to sym- 
metry subspaces [15]. More importantly, there arc several properties of real systems that 
require explicit inclusion of symmetries and they are defined by a variety of Lie algebras. 
This led to studies on EGEs with symmetries such as spin [16-20], spin-isospin SU{4) [21], 
J-symmetry [22] and many others (see for example [23, 24]). In the present paper, we con- 
sider parity symmetry in EE as there are several nuclear structure quantities that require 
explicit inclusion of parity. Some of these are as follows. 
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Parity ratios of nuclear level densities is an important ingredient in nuclear astrophysical 
applications. Recently, a method based on non-interacting Fermi-gas model for proton- 
neutron systems has been developed and the parity (tt) ratios as a function of excitation 
energy in large number of nuclei of astrophysical interest have been tabulated [25]. The 
method is based on the assumption that the probability to occupy s out of given single 
particle (sp) states follow Poisson distribution in the dilute limit {m « N, N ^ oo where 
m is the number of particles) . Then the ratio of the partition functions for the -|-ve and — ve 
parity states is given by the simple formula = tanh/, where / is average number 

of particles in the -|-ve parity states. Starting with this, an iterative method is developed 
with inputs from the Fermi-Dirac distribution for occupancies including pairing effects and 
the Fermi-gas form for the total level density. In the examples studied in [25], parity ratios 
are found to equilibrate only around 5 — 10 MeV excitation energy. However, ab-initio 
interacting particle theory for parity ratios is not yet available. 

A closely related question is about the form of the density of states defined over spaces 
with fixed-TT. In general, fixed-7r density of states can be written as a sum of appropriate 
partial densities. In the situation that the form of the partial densities is determined by a 
few parameters (as it is with a Gaussian or a Gaussian with one or two corrections), it is 
possible to derive a theory for these parameters and using these, one can construct fixed- tt 
density of states and calculate parity ratios. Such a theory with interactions in general 
follows from random matrix theory [15]. 

In addition to the questions related to fixed- tt density of states and parity ratios, there 
is also the important recognition in the past few years that random interactions generate 
regular structures [22, 26, 27]. It was shown in [28] that shell model for even-even nuclei 
gives preponderance of -|-ve parity ground states. A parameter-free EGOE with parity has 
been defined and analyzed recently [29] to address the question of 'preponderance of ground 
states with positive parity' for systems with even number of fermions. They show that in the 
dilute limit, -|-ve parity ground states appear with only 50% probability. Thus, a random 
matrix theory describing shell model results is not yet available. 

With the success of the embedded random matrix ensembles (EE) [9, 14], one can argue 
that the EE generated by parity preserving random interaction may provide generic results 
for the three nuclear structure quantities mentioned above. For nuclei, the GOE versions of 
EE are relevant. Then, with a random (modeled by GOE) two-body interaction preserving 
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parity in the presence of a mean-field, we have embedded Gaussian orthogonal ensemble 
of one plus two-body interactions with parity [hereafter called EGOE(l-|-2)-7r]. This model 
contains two mixing parameters and a gap between the -|-ve and — ve parity sp states and it 
goes much beyond the simpler model considered in [29]. In the random matrix model used 
in the present paper, proton-neutron degrees of freedom and angular momentum (J) are not 
considered. Let us add that in the present paper for the first time a random matrix theory 
for parity ratios is attempted. Now we will give a preview. 

Section II gives the definition of EGOE(l-|-2)-7r and a method for its construction. Prom 
the results known for EE for spinless fermion (boson) systems, for fermions (bosons) with 
spin and from shell model calculations [3, 14, 18], it is expected that the fixed- tt state den- 
sities (more appropriately partial densities) approach Gaussian form in general. Therefore, 
exact propagation formulas for fixed- tt energy centroids and spectral variances are derived 
and the results are given in Section III. Used here is the group theoretical formulation devel- 
oped by Chang et al [30]. Similarly in Appendix B, given are the formulas for the ensemble 
averaged skewness 71 (mi, 7712) and excess 72(^1,7712) parameters for fixed-(777i, 7^2) partial 
densities (with mi fermions distributed in A^_|_ number of +ve parity sp levels and similarly 
1712 fermions in number of — ve parity sp levels) and used is the binary correlation approx- 
imation method described in [31-34]. These will provide corrections to the Gaussian state 
densities. In Section IV, presented are the numerical results for (i) fixed- tt state densities, 
(ii) parity ratios of state densities and (iii) probability for -|-ve parity ground states. Finally, 
Section V gives conclusions and future outlook. 

II. EGOE(l-|-2)-7r ENSEMBLE 

Given A^_|_ number of positive parity sp states and similarly N_ number of negative parity 
sp states, let us assume, for simplicity, that the -|-ve and — ve parity states are degenerate and 
separated by energy A (see Fig. 1). This defines the one-body part h{l) of the Hamiltonian 
H with N = + N_ sp states. The matrix for the two-body part V{2) of H [we assume H 
is (l-|-2)-body] will be a 3 x 3 block matrix in two particle spaces as there are three possible 
ways to generate two particle states with definite parity: (i) both fermions in -|-ve parity 
states; (ii) both fermions in — ve parity states; (iii) one fermion in -|-ve and other fermion 
in — ve parity states. They will give the matrices A, B and C respectively in Fig. 1. For 



4 



parity preserving interactions only the states (i) and (ii) will be mixed and mixing matrix is 
D in Fig. 1. Note that the matrices A, B and C are symmetric square matrices while D is 
in general a rectangular mixing matrix. Consider N sp states arranged such that the states 
1 to A^_|_ have +ve parity and states A^_|_ + 1 to iV have — ve parity. Then the operator form 
of H preserving parity is, 

H = h{l) + V{2) ; 
HI) - Y.et'^n^^ + f: et'nt'; = , = A , 

i=l i=N++l 
N+ 

V{2) = ^ {i^k i^e\V \ui Uj) al al aj ai 

i,j,k,l = 1 
{i <j,k< e) 



N 



(1) 

{i' < f, k' < f ) 
+ ^ ^ (i^fc" I V I Vii' i/j") al„ a\„ a^// Oj// 

i",k"=l j"/"=N++l 

N+ N 

+ 2Z ^ {^p^q\V \ Ps) dp as an + h.c] . 

P,Q = l R,S = N+ + 1 
(P <Q) {R< S) 

In Eq. (1), z/j's are sp states with i = 1,2,. . . ,N (the first A^_|_ states are +ve parity and 
remaining — ve parity). Similarly, {. . . \ V \ . . .) are the two-particle matrix elements, fii are 
number operators and a] and are creation and annihilation operators respectively. Note 
that the four terms in the RHS of the expression for V{2) in Eq. (1) correspond respectively 
to the matrices A, B, C and D shown in Fig. 1. 

Many particle states for m fermions in the N sp states can be obtained by distributing 
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FIG. 1. Parity preserving one plus two-body H with a sp spectrum defining h{l) along with a 
schematic form of the V{2) matrix. Dimension of the matrices A, B and C are A^+(A^+ — l)/2, 
N^{N- — l)/2 and iV_|_iV_ respectively. Note that D-'^ is the transpose of the matrix D. See text 
for details. 

mi fermions in the +vc parity sp states (A^_|_ in number) and similarly, 1712 fermions in the 
— ve parity sp states {N^ in number) with m = mi + m2- Let us denote each distribution 
of mi fermions in sp states by mi and similarly, m2 for m2 fermions in sp states. 
Many particle basis defined by (mi, 1112) with m2 even will form the basis for +ve parity 
states and similarly, with m^. odd for — ve parity states. In the (mi, m2) basis with m2 even 
(or odd) , the H matrix construction reduces to the matrix construction for spinless fermion 
systems. The method of construction for spinless fermion systems is well known [14] and 
therefore it is easy to construct the many particle H matrices in +ve and — ve parity spaces. 
The matrix dimensions 0?+ for +ve parity and d- for — ve parity spaces are given by, 

'"""in. E nr- 

mij\m.2j ^ VmiyVm2 



(2) 



Vmi/ \m2j ^-^ 

m\,m2 (m2 even) ^ ' ^ ' mi,m2 (m2 odd) 

Some examples for the dimensions d+ and d- are given in Table I. 

The EGOE(l+2)-7r ensemble is defined by choosing the matrices A, B and C to be 
independent GOE's with matrix elements variances f^, vl and respectively. Similarly the 
matrix elements of the mixing D matrix are chosen to be independent (independent of A, B 
and C matrix elements) zero centered Gaussian variables with variance v^. Without loss of 
generality we choose A = 1 so that all the v's are in A units. This general EGOE(l+2)-7r 
model will have too many parameters (v^, v'l, v^, vj, N^, m) and therefore it is necessary 



TABLE I. Hamiltonian matrix dimensions d+ and ri_ for various values of {N^,N-,m). 
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to reduce the number of parameters. A numerically tractable and physically relevant (as 
discussed ahead) restriction is to choose the matrix elements variances of the diagonal blocks 
A, B and C to be same and then we have the EGOE(l+2)-7r model defined by (A^+, N_,m) 
and the variance parameters (r,Q;) where 

'^g _ ^ ^ ^ ^ 2 ^ _ 2 /o\ 

A2 A2 A2 ' A2 ^ 
Thus EGOE(l+2)-7r we employ is 

A : GOE(0 -.T^) , B: GOE(0 : r^) , C : GOE(0 : t^) , D : GOE(0 : a^) ; 
A, 5, C, D are independent GOE's . 



(4) 



Note that the D matrix is a GOE only in the sense that the matrix elements Dij are all 
independent zero centered Gaussian variables with variance a^. In the limit — )■ cxd and 
a — T, the model defined by Eqs. (1), (3) and (4) reduces to the simpler model analyzed in 
[29]. 

Before proceeding further, it is useful to mention that we are considering in this paper 
spinless fermion systems (with parity included) just as in the previous investigation [29] . It 



is possible to extend the ensemble to nucleons in shell model j-orbits (including both +ve 
and — ve orbits such as for example sd and fp) and construct the ensemble in many-nucleon 
spaces with a given J'^ or J'^T using a shell model code. However, such an attempt has not 
been made, just as in [29], as our focus is on parity. Also the ensemble for spinless systems 
will give the essential features due to parity and these can be used in later explorations using 
ensembles with J'^ or J'^T which are more complicated numerically and more importantly 
from analytical point of view [3, 22]. In fact, due to the severe problems associated with 
analytical tractabihty, a variety of EGOE are being analyzed since 1995; see [9, 14, 22] for 
reviews. At this point it is also useful to mention that EG0E(l+2)'s are also called TBRE in 
literature; see Section 5.7 in [9] for clarifications on this nomenclature. As Brody et al state 
[3]: The most severe mathematical difficulties with TBRE are due to angular momentum 
constraints . . . Another type of ensemble, . . . much closer to being mathematical tractable 
abandons the J restrictions entirely . . . an embedded GOE, or EGOE for short. 

Starting with the EGOE(l+2)-7r ensemble defined by Eqs. (1), (3) and (4), we have 
numerically constructed 100 (in some examples 200) members of the ensemble in many- 
particle +ve and — ve parity spaces with dimensions dj^ and d^ given by Eq. (2) for several 
values of (iV+, iV_, m) and varying the parameters r and a. This means we have considered 
100 realizations of EGOE(l-|-2)-7r random matrices in (A^+, A^_, m) spaces - we use the phrase 
'members' throughout the paper instead of 'realizations' (other names used by some authors 
are 'sets', 'samples' and 'trials') as in all our previous papers. Before discussing the numerical 
calculations, we present the results for the energy centroids, variances and also the shape 
parameters (skcwness and excess) defining the normalized fixed-(mi,m2) partial densities 



where, (. . .) corresponds to average and a and [5 are extra labels required to specify com- 
pletely the states with a given (mi, m2) and E respectively. Later we use the symbol ((...)) 
that denotes the trace. These will allow us to understand some of the numerical results. 
Let us add that the fixed-7r eigenvalue densities I±{E) are sum of the appropriate partial 
densities as given by Eq. (16) ahead. Note that the densities I±{E) are normalized to d±. 




(5) 



mi, m2, a) = X] Cg^^^'™"'" \E, (5) , 
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III. ENERGY CENTROIDS, VARIANCES, SKEWNESS AND EXCESS PARAM- 
ETERS FOR FIXED- (mi, ma) PARTIAL DENSITIES 



Let us call the set of +ve parity sp states as unitary orbit ^1 and similarly the set of — ve 
parity sp states as unitary orbit #2; see [15] for unitary orbits notation and significance. 
For convenience, from now on, we denote the sp states by the roman letters . . .) and 
unitary orbits by greek letters {a, (3, . . .). Note that a — 1 corresponds to the +ve parity 
unitary orbit and a — 2 corresponds to the — ve parity unitary orbit (with this notation, 
Ni = and A^2 = -^-)- The sp states that belong to a unitary orbit a are denoted as 
iai jai ■ ■ ■■ Propagation formulas for the energy ccntroids and variances of the partial densities 
pmi,m2^^^ follow from the unitary decomposition of V{2) with respect to the sub-algebra 
U{N^)®U{NS) contained in U{N). Note that (mi, ma) label the irreducible representations 
(irreps) of U {N^)®U {NJ) and they all belong to the U (N) irreps labeled by m. The (mi, ma) 
are often called unitary configurations [15]. With respect to U{N^) © U{N_), the operator 
V{2) decomposes into three parts V{2) -> + V^'^^ + I^^l The V^^^ generates the energy 
centroids (y^"^^'"^^^ corresponds to the 'algebraic' mean-field generated by V and V^'^^ 
is the remaining irreducible two-body part. Extending the unitary decomposition for the 
situation with a single orbit for spinless fermions (given in Appendix A) and also using the 
detailed results in [30], we obtain the following formulas for the F^'^l's. The V^°^ is given by 
(with q; = 1, 2 and ^ = 1, 2) 





i>j 



(6) 
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Then the traceless part V - = yl^l + ^I^] ^here = Vic^jpicip - K/3 and 

(y)ijU — Vijke for all others. Now the V^^^ part is 



(7) 



Finally, the part is as follows, 



. . = i/. . . . 



•^fcai^fcci^ ~ ^kcipkaJn AT _ 9A . ' T= ' 



(8) 



= fo^ all others . 



Given Eqs. (6), (7) and (8), by intuition and using Eq. (A3), it is possible to write the 
propagation formulas for the energy ccntroids and variances of p"^'^'"^'^{E). Note that these 
are essentially traces of H and over the space defined by the two-orbit configurations 
(mi, m2); see Eqs. (9) and (10) ahead. A direct approach to write the propagation formulas 
for centroids and variances for a multi-orbit configuration was given in detail first by French 
and Ratcliff [35] . The formula for the variance given in [35] is cumbersome and it is realized 
later [30] that they can be made compact by applying group theory (see also [14, 15, 36]). 
We have adopted the group theoretical approach for the two-orbit averages and obtained 
formulas. Propagation formula for the fixed- (mi, 777.2) energy centroids is, 

^e(mi, m2) = {Hr'""^ = m2 A + '^j'^^-^-^^ V^^ . (9) 
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First term in Eq. (9) is generated by h{l) and it is simple because of the choice of the sp 
energies as shown in Fig. 1. Propagation formula for fixed- (mi, 7712) variances is, 



mi,m2l2 



(10) 



The 'prime' over summations in Eq. (10) implies that the summations are not free sums. 
Note that (a, /3, 7, 5) take values (1, 1, 1, 1), (2, 2, 2, 2), (1, 2, 1, 2), (1, 1, 2, 2) and (2, 2, 1, 1). 
Similarly, in the sum over (zq,, j^), i < j' if a = ^ and otherwise the sum is over all i and j. 
Similarly, for {k^,^s). Using Ec{mi,m2) and (T^(mi,m2), the fixed-parity energy centroids 
and spectral variances [they define I±{E)] can be obtained as follows, 

1 ' 

Ec{m,±) = {H)""'^ = — d{mi,m2)Ec{mi,m2) , 

"± ^ 

mi,m2 

^jj2^m,±^l [a2(r„,,r„2) + £;c'K,^2)] ■ 

a± ^-^ 

The 'prime' over summations in Eq. (11) implies that m2 is even (odd) for -|-ve(— ve) parity. 

It should be pointed out that the formulas given by Eqs. (9), (10) and (11) are compact 
and easy to understand compared to Eqs. (10)-(14) of [29] and also those that follow 
from Eqs. (129) and (133) of [35] where unitary decomposition is not employed. We have 
verified Eqs. (9) and (10) by explicit construction of the H matrices in many examples. In 
principle, it is possible to obtain a formula for the ensemble averaged variances using Eq. 
(10); the ensemble averaged centroids derive only from h{l). Simple asymptotic formulas 
for ensemble averaged variances follow by neglecting the 5-functions that appear in Eq. (10) 
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and replacing V^f.^ by and appropriately. Then the final formula for the ensemble 
averaged fixed-(mi, 771.2) variances is, 

2 

^ m« {N^ - ma) 



(7^ (mi, 7712) ~ m 



a=l 
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nil 
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Here, mi = Ni—mi and m2 = N2 — m,2. The 'overline' in Eq. (12) denotes ensemble average. 
In Table II, we compare the results obtained from Eq. (12) with those obtained for various 
100 member ensembles using Eq. (10) and the agreements are quite good. Therefore, in 
many practical applications, one can use Eq. (12). 

The skewness and excess parameters 71 and 72 give information about the shape of the 
partial densities and they are close to zero implies Gaussian form. Formulas for the Mr, 
r = 3, 4 for a given one plus two-body Hamiltonian defined by Eq. (1) follow from the results 
given in [36-41] many years back. However, these formulas contain very large number of 
complicated terms (in particular for M4) and carrying out analytically ensemble averaging 
is proved to be impractical (we are not aware if anyone was successful in the past). Some 
idea of the difficulty in carrying out simplifications can be seen from the attempt in [42] . An 
alternative is to program the exact formulas and evaluate the moments numerically for each 
member of EGOE(l+2)-7r by considering say 500 members in two particle spaces. As pointed 
out by Teran and Johnson [43] in their most recent attempt, these calculations for the 4th 
moments are time consuming if not impractical. All the problems with the exact formulas 
have been emphasized in [15]. Because of these (in future with much faster computers it may 
be possible to use the exact formulas), we have adopted the binary correlation approximation, 
first used by Mon and French [31, 32] and later by French et al [33, 34], which is good in 
the dilute limit: mi, Ni,m2, N2 00, m/Ni — > and m/A^2 0, where m is mi or m2, 
for deriving formulas for the ensemble averaged M3 and M4. The final formulas are given 
in Appendix B and details of the derivations will be reported elsewhere [44] . The following 
results are inferred from the results in Appendix B. 

It is seen from Eq. (B9), 7i(mi,m2) will be non-zero only when a 7^ and the r 
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TABLE II. Ensemble averaged fixed-(mi, 7712) widths a{mi,m2) and the total spectral width at 
for different (r, a) values. For each (r, a), the cr(mi, 777-2) are given in the table and they are 
obtained using the exact propagation formula Eq. (10) for each member of the ensemble. In all the 
calculations, ensembles with 100 members are employed. Numbers in the bracket are obtained by 
using the asymptotic formula given in Eq. (12). Last row for each (A?+, 7V_) gives the corresponding 
at values. All the results are given for 6 particle systems and the dimensions (i(mi, 777-2) are also 
given in the table. See text for details. 
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dependence is weak. Also, it is seen that for — N^, 71 (mi, 777,2) = —71 (^2, mi). Similarly, 
Eq. (BIO) shows that for — N^, 72(7771, 7772) = 72(m2, 7?7i). In the dilute limit, with some 
approximations as discussed after Eq. (BIO), the expression for 72(7712,^1) is given by Eq. 
(Bll). This shows that, for a << r or t << a, the Ci and C2 in Eq. (Bll) will be negligible 
and then, 72 ~ — 4/m for rrii = 77^2 =777/2 and Ni = N2 = N. This is same as the result for 
spinless fermion EG0E(2) [31, 32] and shows that for a range of (r, a) values, p^^''^^{E) will 
be close to Gaussian. Moreover, to the extent that Eq. (Bll) applies, the density p^'^''^'^(E) 
is a convolution of the densities generated by ^(2) and -0(2) operators. Let us add that the 
binary correlation results presented in Appendix B, with further extensions, will be useful 
in the study of partitioned EGOE discussed in [14, 45]. 



IV. RESULTS AND DISCUSSION 

In order to proceed with the numerical calculations, we need to have some idea of the 
range of the parameters {T,a,m/N^, N^/N^). Towards this end, we have used realistic 
nuclear effective interactions in sdfp [46] and fpgg/2 [47] spaces and calculated the variances 
vl, vl, v1, v1 for these interactions. Note that it is easy to identify the matrices A, B, 
C and D given the interaction matrix elements {{3132) JT \ V \ {3z3a)JT). To calculate the 
mean-squared matrix elements I'^'s, we put the diagonal two-particle matrix elements to 
be zero and use the weight factor (2 J -|- 1)(2T -|- 1). Assuming that A = 3 MeV and 5 
MeV (these are reasonable values for A = 20 — 80 nuclei), we obtain r ~ 0.09 — 0.24 and 
a ~ (0.9 — 1.3) X r. These deduced values of a and r clearly point out that one has to 
go beyond the highly restricted ensemble employed in [29] and it is necessary to consider 
the more general EGOE(l-|-2)-7r defined in Section II. Similarly, for sdfp and fpgg/2 spaces 
Nj^/N_ ~ 0.5 — 2.0. Finally, for nuclei with 777 number of valence nucleons (particles or holes) 
where sdfp or fpgg/2 spaces are appropriate, usually 777 ^ or A?^., whichever is lower. 
Given these, we have selected the following examples: (A^+, A^_, 777) = (8,8,4), (8,8,5), 
(10, 6, 4), (10, 6, 5), (6, 10, 4), (6, 10, 5), (8, 8, 6), (6, 6, 6), (7, 7, 7) and (7, 7, 6). To go beyond 
the matrix dimensions ~ 5000 with 100 members is not feasible at present with the HPC 
cluster that is used for all the calculations. Most of the discussion in this paper is restricted 
to A" = A^+ -I- A^_ = 16 and 777 << A" as in this dilute limit it is possible to understand the 
ensemble results better. Following the nuclear examples mentioned above, we have chosen 
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T = 0.05, 0.1, 0.2, 0.3 and cc/r = 0.5, 1.0, 1.5. We will make some comments on the results 
for other (r, a) values at appropriate places. 

Now we will present the results for (i) the form of the +ve and — ve parity state densities 

I+{E) and I~-{E) respectively, (ii) the parity ratios I^{E)/ I^{E) vs E where E is the exci- 
tation energy of the system and (iii) the probability for -|-ve parity ground states generated 
by the EGOE(l+2)-7r ensemble. 

A. Gaussian form for fixed-7r state densities 

Using the method discussed in Sec. II, we have numerically constructed in -|-ve and — ve 
parity spaces EGOE(l+2)-7r ensembles of random matrices consisting of 100 Hamiltonian 
matrices in large number of examples, i.e. for (A^+, N^, m) and (r, a) parameters mentioned 
above. Diagonalizing these matrices, ensemble averaged eigenvalue (state) densities, 

l^) = mH-E))t, (13) 

arc constructed. From now on, wc drop the 'ovcrline' symbol when there is no confusion. 
Results are shown for (iV+,iV_,m) = (8,8,4), (8,8,5), (10,6,5) and (6,10,5) for several 
values of (r, a) in Figs. 2, 3 and 4. To construct the fixed-parity eigenvalue densities, 
we first make the centroids Ec{m, ±) of all the members of the ensemble to be zero and 
variances (7^(m,±) to be unity, i.e. for each member we have the standardized eigenvalues 
E — [E — Ec{m, ±)]/(T(m, ±). Then, combining all the E and using a bin-size AE — 0.2, 
histograms for I±{E) are generated. It is seen that the state densities are multimodal for 
small r values and for r > 0.1, they are unimodal and close to a Gaussian. Note that in our 
examples, a — (0.5 — 1.5) x r. 

For V{2) — 0, the eigenvalue densities will be a sum of spikes at 0, 2A, 4A, ... for -|-ve 
parity densities and similarly at A, 3A, 5A, ... for — ve parity densities. As we switch on 
y(2), the spikes will spread due to the matrices A, B and C in Fig. 1 and mix due to the 
matrix D. The variance cr^ (mi, 777.2) can be written as, 

(7^(r77i, 7772) = (T^(777i, 7772 — > '>Tli,m2) + (T^(777i, 7772 —>■ 777 1 ± 2, 7772 ^ 2) ■ (14) 

The internal variance (7^(7771,7772 — > 7771,7772) is due to A, B and C matrices and it receives 
contribution only from the r parameter. Similarly, the external variance (7^(7771,7772 — > 
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FIG. 2. (Color online) Positive and negative parity state densities for various (r, a) values for 
(A^+, A^_,m) = (8,8,4) system. See text for details. 
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(N,,N_,m)=(8,8,5) 
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FIG. 3. (Color online) Positive and negative parity state densities for various (r, a) values for 
{Nj^,N^,m) = (8,8,5) system. Histograms are numerical ensemble results. The dashed (red) 
curve corresponds to Gaussian form for p"^'^''^^(E) in Eq. (16) and similarly, solid (green) curve 
corresponds to Edgeworth corrected Gaussian form with 71 (mi, 777-2) and 72(^,1, 7712) obtained using 
the results in Appendix B. See text for details. 
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FIG. 4. (Color online) Positive and negative parity state densities for various (r, a) values for 
{Nj^,N^,m) = (10,6,5) and (6,10,5) systems. Histograms are numerical ensemble results. The 
dashed (red) curve corresponds to Gaussian form for p"^'^''^^(E) in Eq. (16) and similarly, solid 
(green) curve corresponds to Edgeworth corrected Gaussian form with 7i(mi,m2) and 72(^7^l,m2) 
obtained using the results in Appendix B. See text for details. 
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FIG. 5. (Color online) Positive and negative parity state densities for some small values of (r, a). 
The {N-^-,N-,m) values are given in the figures. See text for details. 



mi ±2, 777,2 =f2) is due to the matrix D and it receives contribution only from the a parameter. 
When we switch on V{2), as the ensemble averaged centroids generated by V{2) will be zero, 
the positions of the spikes will be largely unaltered. However, they will start spreading and 
mixing as r and a increase. Therefore, the density will be multimodal with the modes well 
separated for very small (r, a) values. Some examples for this are shown in Fig. 5. As r and 
a start increasing from zero, the spikes spread and will start overlapping for cr{mi,'m2) ^ A. 
This is the situation with r = 0.05 shown in Figs. 2, 3 and 4. However, as r increases (with 
a ~ r), the densities start becoming unimodal as seen from the r = 0.1 and 0.2 examples. 
Also, the m dependence is not strong as seen from the Figs. 2, 3 and 4. Now we will 
discuss the comparison of the ensemble results with the smoothed densities constructed 
using ^'^(mi, 7772), cr^(mi, 7772), 71(^1,7772) and 72(^1,7772). 
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FIG. 6. (Color online) Positive and negative parity state densities for various (r, a) values for 
{Nj^, N^,m) = (8,8,6), (6,10,6) and (10,6,6) systems. Smoothed curves (solid red lines) are 
obtained using fixed-(mi, 771,2) partial densities. See text for details. 
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As the particle numbers in the examples shown in Figs. 2, 3 and 4 are small, the excess 
parameter 72(771,71) ~ —0.7 to —0.8 (skewness parameter 7i(m, tt) ~ in all our examples). 
Therefore the densities are not very close to a Gaussian form. It has been well established 
that the ensemble averaged eigenvalue density takes Gaussian form in the case of spinless 
fermion (as well as boson) systems and also for the embedded ensembles extending to those 
with good quantum numbers; see [3, 9, 14, 18] and references therein. Thus, it can be 
anticipated that Gaussian form is generic for the state densities or more appropriately, for 
the partial densities p^^''^^(E) generated by EGOE(l+2)-7r for some range of (r, a) values. 
Results for the fixed- tt densities for (A^+, N_,m) = (8, 8, 6), (6, 10, 6) and (10, 6, 6) systems 
arc shown in Fig. 6. The smoothed +ve and — ve parity densities arc a sum of the partial 
densities p"i'"2(^), 

P^^^) -l-Y. dim,,m2)p^^''^'{E) . (15) 



mi, 1712 



Note that the summation in Eq. (15) is over m2 even for -|-ve parity density and similarly 
over 1712 odd for — ve parity density. Here p±{E) as well as p^'^''^'^(E) are normahzed to 
unity. However, in practice, the densities normalized to dimensions are needed and they are 
denoted, as used earher, by I±{E) and I'^'^^'^'^i^E) respectively, 

/ 

I±{E) = d±p±{E) = /-i'-2(E) ; jmum^l^^^ ^ d{mi , 7112) p""' '"^^ (E) • (16) 

mi,m2 

We employ the Edgeworth (ED) form that includes 71 and 72 corrections to the Gaussian 
partial densities p^^'"^^{E). Then 

and in terms of the standardized variable E, the ED form is given by 

riED{E) = r^g{E) |l + \^-^He^ (e)] + V^^He^ (e) + ^^He, {e) 



1 I E"^ 



(17) 



Here, He are Hermite polynomials: Hez{x) — — 3x, Hei{x) — x'^ — Qx^ + 3 and HeQ{x) — 
x^ — 15x^ + 45a;^ — 15. Using Eqs. (15) and (17) with exact centroids and variances given 
by the propagation formulas in Section HI and the binary correlation results for 71 and 72 
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as given by the formulas in Appendix B, the smoothed +ve and — ve parity state densities 
are constructed. We put r]ED{E) — when riED{E) < 0. It is clearly seen from Fig. 6 
that the sum of partial densities, with the partial densities represented by ED corrected 
Gaussians, describe extremely well the exact fixed- vr densities in these examples. Therefore, 
for the (r, a) values in the range determined by nuclear sdfp and fpg9/2 interactions, i.e. 
T ~ 0.1 — 0.3 and a ~ 0.5r — 2r, the partial densities can be well represented by ED 
corrected Gaussians and total densities are also close to ED corrected Gaussians. Unhke 
Fig. 6, densities in Figs. 2, 3 and 4 show, in many cases, strong departures form Gaussian 
form. Therefore, it is important to test how well Eq. (16) with ED corrected Gaussian 
for p™i'™2(i^) describes the numerical results for I±{E). We show this comparison for all 
the densities in Figs. 3 and 4. It is clearly seen that the agreements with ED corrected 
Gaussians are good in all the cases. Therefore, the large deviations from the Gaussian form 
for I±{E) arise mainly because of the distribution of the centroids [this involves dimensions 
of the (777,1,777,2) configurations] of the partial densities involved. It is possible that the 
agreements in Figs. 3 and 4 may become more perfect if we employ, for the partial densities, 
some non-canonical forms defined by the first four moments as given for example in [48, 49]. 
However, as these forms are not derived using any random matrix ensemble, we haven't used 
these for the partial densities in our present investigation. In conclusion, for the physically 
relevant range of (r, a) values, the propagation formulas for centroids and variances given 
by Eqs. (9) and (10) or alternatively with Ec{mi,m2) — and Eq. (12) along with the 
EGOE(l-|-2)-7r ensemble averaged 71 (mi, 7712) and 72(7771,7772) formulas (obtained using the 
binary correlation approximation as given in Appendix B) can be used to construct fixed- tt 
state densities for larger {N^, N_,m) systems. 

B. Parity ratios for state densities 

As stated in the beginning, parity ratio of state densities at a given excitation energy (E) 
is a quantity of considerable interest in nuclear structure. For the systems shown in Figs. 2, 
3 and 4 and also for many other systems, we have studied the parity ratios and the results are 
shown in Figs. 7-10. As the parity ratios need to be calculated at a given value of excitation 
energy E, we measure the eigenvalues in both -|-ve and — ve parity spaces with respect to 
the absolute ground state energy Egg of the N — A^_,_ -|- N_ system. Thus, Egg is defined by 
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taking all the +ve and — ve parity eigenvalues of all members of the ensemble and choosing 
the lowest of all these. The ground state energy can also be determined by averaging the +ve 
and — ve parity ground state energies over the ensemble and then the ground state energy is 
minimum of the two. It is seen that the results for parity ratios are essentially independent 
of the choice of Egg and thus we employ absolute ground state energy in our calculations. 
We use the ensemble averaged total (+ve and — ve eigenvalues combined) spectrum width 
at of the system for scahng. The total widths at can be calculated also by using Ec{mi,m2) 
and (T^ (mi, 7712). Examples for at are shown in Table II and they are in good agreement 
with the results obtained using the simple formula given by Eq. (12). We use the variable 
E = {E — Egs)/at for calculating parity ratios. Starting with Egg and using a bin-size of 
AE = 0.2, wc have calculated the number of states I+{E) with +ve parity and also the 
number of states I-{E) with — ve parity in a given bin and then the ratio I-{E)/I^{E) is 
the parity ratio. Note that the results in Figs. 7-10 are shown ior E — — 3 as the spectrum 
span is ~ 5.5(7^. To go beyond the middle of the spectrum, for real nuclei, one has to include 
more sp levels (also a finer splitting of the -|-ve and — ve parity levels may be needed) and 
therefore, and change. Continuing with this, one obtains the Bethe form for nuclear 
level densities [15]. 

General observations from Figs. 7-10 are as follows, (i) The parity ratio I-{E)/I+{E) 
will be zero up to an energy Eq. (ii) Then, it starts increasing and becomes larger than 
unity at an energy Em- (iii) From here on, the parity ratio decreases and saturates quickly 
to unity from an energy Ei. In these examples, Eq ^ 0.4, E„i ~ 1 and Ei ~ 1.5. It 
is seen that the curves shift towards left as r increases. Also the position of the peak 
shifts to much larger value of E^ and equilibration gets delayed as a increases for a fixed 
T value. Therefore for larger r, the energies (£^0, Em, Ei) are smaller compared to those 
for a smaller r. The three transition energies also depend on N_,m). We have also 
verified, as shown in Fig. 9, that the general structure of the parity ratios will remain same 
even when we change A — )> —A (i.e. — ve parity sp states below the +ve parity sp states). 
For {N_^_, N_,m) = (8,8,4) system, results for A = 1 are given in Fig. 7 and they are 
almost same as the results with A = — 1 given in Fig. 9. The general structures (i)-(iii) are 
clearly seen in the numerical examples shown in [25] where a method based on the Fermi-gas 
model has been employed. If 0"^ ~ 6 — 8 MeV, equilibration in parities is expected around 
£■ ~ 8 — 10 MeV and this is clearly seen in the examples in [25] . It is also seen from Fig. 
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FIG. 7. (Color online) Parity ratios for various (r, a) values for (A^+, N-,m) = (8, 8, 4) and (10, 6, 4) 
systems. See text for details. 



8 that the equilibration is quite poor for very small values of r and therefore comparing 
with the results in [25] , it can be argued that very small values of r are ruled out for nuclei. 
Hence, it is plausible to conclude that generic results for parity ratios can be derived using 
EGOE(H-2)-7r with reasonably large (r, a) values. Let us add that the interpretations in 
[25] are based on the occupancies of the sp orbits while in the present work, they are in 
terms of r and a parameters. 

Using the smoothed I^{E), constructed as discussed in Section IV A, smoothed forms 
for parity ratios are calculated as follows. Starting with the absolute ground state energy 
Egs and using a bin-size of AE = 0.2, +ve and — ve parity densities in a given energy bin 
are obtained and their ratio is the parity ratio at a given E. We have chosen the examples 
where /_|_ and /_ are close to Gaussians. It is seen from Fig. 10 that the agreement with 
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FIG. 8. (Color online) Parity ratios for various (t, a) values for (A^+,A^_,m) = (6,10,4) and 
(6, 10, 5) systems. See text for details. 
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A=-l 

(N^,N_,m)=(5,1 0,4) (N^,N_,m)=(8,8,4) 



t=0.1,oc=0.5t 



T=0.1,oc=1.5x 



T=0.2,a=0.5T 



T=0.2,oc=1.5'C 



x=0.1,oc=0.5t 



2 



1 - 



T=0.1,a,= 1.5T 



T=0.2,a=0.5T 



1 - 



T=0.2,a=1.5T 



t 

1 2 3 1 2 3 

(E-E^3)/a, 

FIG. 9. (Color online) Parity ratios for some values of (r, a) with A = —1 for {Nj^., N-,m) 
(5, 10, 4) and (8, 8, 4) systems. See text for details. 



exact results is good for E ^ 0.5. However, for smaller E, to obtain a good agreement one 
should have a better prescription for determining the tail part of the p^^'^'^{E) distributions. 
Developing the theory for this is beyond the scope of the present paper as this requires more 
complete analytical treatment of the ensemble. 



C. Probability for +ve parity ground states 

Papenbrock and Weidenmiiller used the r — )■ oo, a = r limit of EGOE(H-2)-7r for several 
{N^, N_,m) systems to study the probability for +ve parity ground states over the 
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1=0.2 a=0.5x 



(N^,N_,m)=(1 0,6,6) 
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T=0.2 
a=1.5T 




(N^,N_,m)=(1 0,6,6) 

T=0.1 a=1.5'c 




(N^,N_,m)=(6,10,6) 
t;=0.1 a=0.5x 



1 



1 



(E-E,3)/a, 

FIG. 10. (Color online) Parity ratios for various (r, a) values and for various {Nj^,N^,m) systems. 
Filled squares (brown) are obtained using fixed-(mi, 7712) partial densities. See text for details. 

ensemble [29]. As stated before, this exercise was motivated by shell model results with 
random interaction giving preponderance of +ve parity ground states [28]. The numerical 
calculations in [29] showed considerable variation (18 — 84%) in _R+. In addition, they gave a 
plausible proof that in the dilute limit [m << (A^+, A^^)], R+ will approach 50%. Combining 
these, they argued that the observed preponderance of +ve parity ground states could be 
a finite size (finite A^+, A^_, m) effect. For the extended EGOE(l+2)-7r considered in the 
present work, where the r — )■ 00 and a = t restriction is relaxed, as we will discuss now, i?+ 
can reach 100%. 

For EGOE(H-2)-7r with r ~ 0, clearly one will get = 100% (for even m and m << 
A^+,iV_) and therefore it is of interest to study R+ variation with (r, a). We have carried 
out calculations using a 200 member ensemble for {N^, N_,m) = (6,6,6) and 100 member 
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a=1.5T 



^ (N^,N.,m)=(6,10,5) 
• (N^,N_,m)=(6,10,4) 
^ (N,,N.,m)=(8,8,5) 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



FIG. 11. (Color online) Probability {R+) for +ve parity ground states for various (r, a) values and 
for various {N+,N-,m) systems. See text for details. 

ensembles for (8, 8, 5), (6, 6, 6), (6, 10, 4) and (6, 10, 5) systems. In these calculations, we use 
a = T and 1.5r. The results are shown in Fig. 11. For a = t, the results are as follows. 
For T ^ 0.04, we have ~ 100% and then starts decreasing with some fluctuations 
between r = 0.1 and 0.2. The origin of these fluctuations is not clear. As r > 1 is not 
realistic, we have restricted the i?+ calculations to r < 1. We see from the figure that 
EGOE(l+2)-7r generates R+ ^ 50% for r < 0.3 independent of {N+,N_,m). Also, R+ 



28 



decreases much faster with r and reaches ~ 30% for r = 0.5 for N^,m.) — (6,6,6). 
For m < (A^+, A^_), the decrease in is slower. If we increase a, from the structure of 
the two-particle H matrix in Fig. 1, we can easily infer that the width of the lowest -|-ve 
parity {1711,1712) unitary configuration becomes much larger compared to the lowest — ve 
parity unitary configuration (see Table II for examples). Therefore, with increasing a we 
expect R+ to increase and this is clearly seen in Fig. 11. Thus o; ^ r is required for i?+ to 
be large. A quantitative description of R+ requires the construction of -|-ve and — ve parity 
state densities more accurately in the tail region and the theory for this is not yet available. 



V. CONCLUSIONS AND FUTURE OUTLOOK 



In the present work, we have introduced a generalized EGOE(l-|-2)-7r ensemble for iden- 
tical fermions and its construction follows from EG0E(l-|-2) for spinless fermion systems. 
Using this generalized EE, we have not only studied as it was done by Papenbrock and 
Weidenmiiller [29] using a simpler two-body ensemble with parity, but also studied the form 
of fixed-TT state densities and parity ratios which are important nuclear structure quantities. 
Numerical examples (see Figs. 2-4 and 6), with the range of the various parameters in the 
model fixed using realistic nuclear effective interactions, are used to show that the fixed- tt 
state densities in finite dimensional spaces are of Gaussian form for sufficiently large val- 
ues of the mixing parameters (r, a). The random matrix model also captures the essential 
features of parity ratios as seen in the method based on non-interacting Fermi-gas model 
reported in [25]. We also found preponderance of +ve parity ground states for r ^ 0.5 and 
a ~ 1.5r. In addition, for constructing fixed- tt Gaussian densities we have derived an easy 
to understand propagation formula [see Eq. (10)] for the spectral variances of the partial 
densities p^^''^'^(E) that generate /+ and Similarly, for calculating the corrections to the 
Gaussian forms, formulas for skewness 71 and excess 72 of the partial densities p^'^''^^(E) are 
derived using the binary correlation approximation (sec Appendix B for the formulas). The 
smoothed densities constructed using Edgeworth corrected Gaussians are shown to describe 
the numerical results for I±{E) [for (r, a) values in the range defined by nuclear sdfp and 
fP99/2 interactions - see beginning of Sec. IV] and also the parity ratios at energies away 
from the ground state. Numerical results presented for parity ratios at lower energies show 
that a better theory for the tails of the partial densities is needed (see Figs. 7-10). Thus, 
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the results in the present paper represent considerable progress in analyzing EGOE(l+2)-7r 
ensemble going much beyond the analysis presented in [29]. 

Results in the present work are largely numerical and they clearly show that developing a 
complete analytical theory, going beyond the results presented in Sections III and Appendix 
B, for EGOE(l+2)-7r is important. In future, it is important to investigate EGOE(l+2)-7r 
for proton-neutron systems and then we will have four unitary orbits (two for protons and 
two for neutrons). In addition, by including non-degenerate +ve and — vc parity sp states the 
model could be applied to nuclei for predicting parity ratios. This extended EGOE(l+2)-7r 
model with protons and neutrons occupying different sp states will be generated by a 10x10 
block matrix for V{2) in two-particle spaces. Therefore, parametrization of this ensemble is 
more complex. Analysis of this extended EGOE(l-|-2)-7r is for future. 
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APPENDIX A 



Let us consider a system of m fermions in sp states with a (l-|-2)-body Hamiltonian 
H — h{l) + V{2) where h{l) — J2i V{2) is defined by the two-body matrix elements 

Vijki — {kl I ^(2) I ij). With respect to the U{N) group, the two-body interaction V{2) can 
be separated into scalar (i/ = 0), effective one-body {u — 1) and irreducible two-body {u = 2) 
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parts [14, 15, 30, 36], 

h{h — 1) 



V' 



u=0 



2 

n-1 
N-2 



kikj 



^ijkl 1 



Yu=2 = Y _ Yv=o _ Y^=^ ^ 

^iFfc' = ^i.-^^ -{N- 2)-\^,k for ^ , 
^ijw^ — ^ijki for ^-ll other cases . 
Similar to Eq. (Al), the h{l) operator will have = 0, 1 parts, 

i 

h''=^ = ^ e]hi , e] = ei-e . 



(Al) 



(A2) 



Then the propagation equations for the m-particle centroids and variances are [14, 15, 30, 36], 



m{N — m) 
N{N - 1) 



(A3) 



m{m - 1){N - m){N - m - 1) //. 2y 

Ar(Ar- i)(Ar-2)(Ar-3) \V ^ 



APPENDIX B 



For the EGOE(l+2)-7r Hamiltonian defined in Eq. (1), we have H = h{l) + V{2) = 
h{l) + X{2) + D{2) with X{2) = ^ © S © C is the direct sum of the spreading matrices A, 
B and C and D(2) = D + D is the off- diagonal mixing matrix as shown in Fig. 1. Here, D 
is the transpose of the matrix D. With the sp energies defining the mean field h{l) as given 
in Eq. (1), the first moment Mi of p'^^'"^'^[E) is trivially. 



M,{m,,m2) = {{h + V)y 



m2 , 



(Bl) 
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TABLE III. Exact results for skewness and excess parameters for fixed-vr eigenvalue densities I±{E) 
compared with the binary correlation results (in the table, called 'Approx'). For exact results, we 
have used the eigenvalues obtained from EGOE(l+2)-7r ensembles with 100 members. The binary 
correlation results are obtained using Eqs. (B1)-(B8) and extension of Eq. (11). See text for 

details^ 

7i(m,7r) 72("i,7i") 



(A^+,iV_,m) (r, a/r) Exact Approx Exact Approx 







TT = + 


TT = — 


TT = + 


TT = — 


TT = + 


TT = — 


TT = + 


TT = — 


(8,8,5) 


(0.05,0.5) 


0.15 


-0.15 


0.15 


-0.15 


-0.52 


-0.52 


-0.52 


-0.52 




(0.05, 1.0) 


0.16 


-0.16 


0.16 


-0.16 


-0.50 


-0.50 


-0.50 


-0.50 




(0.05,1.5) 


0.18 


-0.17 


0.18 


-0.18 


-0.46 


-0.46 


-0.46 


-0.46 




(0.2,0.5) 


-0.03 


0.03 


-0.03 


0.03 


-0.71 


-0.71 


-0.71 


-0.71 




(0.2,1.0) 


-0.01 


0.01 


-0.01 


0.01 


-0.73 


-0.73 


-0.74 


-0.74 




(0.2,1.5) 


0.02 


-0.02 


0.02 


-0.02 


-0.72 


-0.72 


-0.73 


-0.73 


(10,6,5) 


(0.05,0.5) 


-0.06 


0.09 


-0.07 


0.09 


-0.26 


-0.76 


-0.26 


-0.75 




(0.05,1.5) 


-0.04 


0.15 


-0.05 


0.15 


-0.01 


-0.86 


-0.01 


-0.86 




(0.2,0.5) 


0.01 


-0.04 


0.01 


-0.04 


-0.73 


-0.69 


-0.73 


-0.69 




(0.2,1.5) 


0.01 


0.02 


0.01 


0.02 


-0.69 


-0.75 


-0.70 


-0.75 


(6,10,5) 


(0.05,0.5) 


-0.09 


0.07 


-0.09 


0.07 


-0.76 


-0.26 


-0.75 


-0.26 




(0.05, 1.5) 


-0.15 


0.05 


-0.15 


0.05 


-0.86 


-0.01 


-0.86 


-0.01 




(0.2,0.5) 


0.04 


-0.01 


0.04 


-0.01 


-0.68 


-0.73 


-0.69 


-0.73 




(0.2,1.5) 


-0.02 


-0.01 


-0.02 


-0.01 


-0.75 


-0.69 


-0.75 


-0.70 



as ^/j,»')™'i'"*2 _ (7722)^ and ((y)'"!'™^ _ g_ gy extending the binary correlation method to 
traces over two-orbit configurations, we have derived formulas for the second, third and 
fourth order traces giving (mi, 7712), r = 2 — 4. It is important to mention that the 
presence of the mixing matrix D makes the derivations lengthy. Therefore, we give only the 
final formulas in this paper and discuss elsewhere the details of the derivations [44]. The 
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second moment M2 is, 



M2(mi, m2) = {{h + Vf)^-'^^ = ^;,2^mi,m. ^ ^^2^-1,^2 

= (m2)^ + ^"(7711, 7712) + V(mi,m2) + T>(mi, 1712) ; 
^ /y2\'»i."^2 2 f fhi + i\ f mA f 7712 + j\ f m2\ 



(B2) 



/ ~\ mi,m2 

Here, for brevity we have defined X{mi,m2) = l^x^)^''"^^ , V{mi,m2) = (DD) and 

^ / ^ \ rni,rn2 

V{mi,m2) = (DD) . Note that, Eq. (B2) gives the binary correlation formula for 
(mi, 777.2) and it reduces to Eq. (12) as expected. Similarly, the third moment M3 is 



^3(7711, 7712) = {{h + y)3)™i'"'2 = (/,3^™i,-2 ^ 2 (y2^™i- 



/ ~\ mi,m2 I ~ \ mi,m2 



+ (XhX)'^'''^' + {^DhDj +(^DhDj 

= (7712)^ + 3(7772) A'(777i, 7772) + (37772 + 2)1^(7771, 7772) + (3772,2 - 2)15(7711 , 7/12) 

(B3) 

The formula for the fourth moment M4 is. 



M4(7ni,7n2) = {{h + V^r""'' = (/i^)-!'-^ ^ 3 ^^2^-1,^2 ^^2^-1,^2 



mi,m2 / ~ \ m\,m2 



, mi,m2 



+2 (hDhD'j + 2 (hDhD'j + (V^4^™i,rn2 
= (7772)^ + 6(7772)^^^(^1, m2) + [6(7772)2 + 8(7772) + 4] P(777i, 7772) 



(B4) 



+ [6(7772)2 - 8(7772) + 4] P(777i, 7772) + 
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The only unknown in Eq. (B4) is (y^)"^^''^'^ and the expression for this is comphcated, 



^\ mi,m2 / ~ \ mi,m2 1 / ^\ 

^y4p,"^2 ^ ^^4^mi,m2 ^ 3 ^^2^mi,m2 ^ ^^^\ + ( DD) [ + ( DXW) 



mi,m2 I / ^\ mi,TO2 



^{DX'^D) +2{XDXD) +2(XDXD) + {{D + D)^ 



= 2[A'(mi, 777.2)]^ + 3[A'(mi,m2)][X>(mi,TO2) + X>(mi, 771.2)] 

+Ti(mi, 7712) + T2 (7771, 7772) + 2X3(7771, 7772) + T4(777i, 7772) . 
where 

Ti(777i,7772) = F {uii, N^, i, t) F {m2, N2, j, u) ] 

1+3=2, t+u=2 



(B5) 



F(7n,7V,^i,fc2) = £(^_ 

s=0 ^ ^ 



777 — s\ fN — m + ki — s\ I m — s\ ( N — m\ (m\ (N + 1 



X 



A^-2s + l fN-s\'^ fk2\~^ 



72("7i, 7772) = 'D{mi, m2)X{mi - 2, 7772 + 2) + V{mi, m2)X{mi + 2, 7772 - 2) , 



73(^1,7772) ^2 ^ 

i+i=2 



mi — i \ ( 777.2 ~J \ _|_ I iTT-i — i \ I 'm'2 — J 



X 



T4(777i,7772) = [1^(7771, 7772)]^ + 1^(7771, 7^2) 



+V{mi, 7772) 2 r»(777i - 2, 7772 + 2) + ^(mi - 2, 7772 + 2) 



+P(777i, 7772) 2 P(777i + 2, 7772 - 2) + ^(^i + 2, 7772 - 2)] + 4 P(777i, 77^2)1^(^1, 7^2) ■ 



(B6) 
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Given the moments Mj.(mi, 777.2) = (^}jr^'^i'i^'2.^ r = 1 — 4, the skewness and excess parameters 
7i and 72 are as follows [50] , 



71(7771,7772) = 



^3(7771, 7772) 

[A;2(777i,7772)]^/^ 



72(7771,7772) = 



A;4 (7771, 7772) 
[A;2 (7771, 7772)] 



2 ' 



(B7) 



where. 



A;2 (7771, 7772) = M2(777i, 7772) - (7771, 7772) , 

/C3 (7771, 7772) = M3(777i, 7772) - 3 M2(777i, 7772) Ml (7771, 7772) + 2 (7771,7772) , 

A;4 (7771, 7772) = M4(777i,7772) - 4 M3(777i,7772) Mi(777i,7772) - 3 M|(777i,7772) 

+ 12 M2(777i,7772) Mf {171^,1712) - 6 M^{mi,m2) . 

After carrying out the simplifications using Eqs. (B1)-(B8), it is easily seen that. 



(B8) 



71(7771,7772) 



V{mi, 7772) - ^^(777i, 7772) 



3/2 



|X>(7771, 7772) + I>(777i, 7772) + ^"(7771, 7772)| 

The expression for 72 is more complex, 

72(7771,7772) = 1^(7771,7772) + 1^(7771,7772) + ^^(7771,7772)! 

X [Ti (7771, 7772) + ^2(7771, 7772) + 2 r3 (777i,7772) + ^4(7771 , 7772) 

I>(777i, 7772) + ©(7771, 7772) [4 - ^"(7771, 7772)] - 2 D(777i, 7772) + '^{mi, 7772 



(B9) 



+ I>(777i,7772) +X'(777i,7772)J [4 - ^"(7771,7772)] - 2 V{mi, 1712) + T>{mi, 1712) "l- 

(BIO) 

With Ti ~ [A" (7771, 7772)]^ + Ci (7771, 7772), T2 = T3 ~ X {mi, 1712) [V {mi, 1712) + V{mi,m2)] and 
T4 ~ 3 [1^(7/71, 7/72) + 2^(7r7i, 7/72)]^ + 6*2(7771, 7772) which are good in the dilute limit (|C1| and 
IC2I will be close to zero), we have 



72(7771,7772) = 



Ci(777i, 7772) + 6*2(7771, 7772) + 4 [1^(7771, 7772) + V{mi,m2)] 

115(7771,7772) + ©(7771,7772) + ^"(7771,7772)1 



(Bll) 



Note that C*i and X depend only on r. Similarly, C*2 and {V,V) depend only on a. The 
{V + V) term in the numerator will contribute to 72(^1,^2) when r = and a is very 
small. The approximation T2 — ^ X{T> + T>) is crucial in obtaining the numerator in 
Eq. (Bll) with no cross-terms involving the a and r parameters. With this, we have to 
be the sum of A;4's coming from X{2) and -0(2) matrices [note that, as mentioned before, 
X{2) = A e S e C and ^(2) ^ D + D]. 
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To test the accuracy of the formulas for Mr given by Eqs. (B1)-(B6), the binary cor- 
relation results for 7i(m,±) and 72(771, ±) are compared with exact results obtained us- 
ing the eigenvalues from EGOE(l+2)-7r ensembles with 100 members for several values of 
(A^+, N_,m) and (r, a) parameters in Table III. Extension of Eq. (11) along with the results 
derived for Mr{mi,m2) will give the binary correlation results for 71(771, ±) and 72(7/7, ±). It 
is clearly seen from the results in the Table that in all the examples considered, the binary 
correlation results are quite close to the exact results. Similar agreements are also seen in 
many other examples which are not shown in the table. 
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